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We discuss alternative homogeneous electron gas systems in which a finite number n of electrons are confined 
to a £)-dimensional sphere. We derive the first few terms of the high-density {rg — > 0, where is the Seitz 
radius) energy expansions for these systems and show that, in the thermodynamic limit (n — >■ cxd), these terms 
become identical to those of J5-dimensional jellium. 
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I. INTRODUCTION 

The Z?-dimensional uniform electron gas (UEG), or 
D-jellium, is the foundation of most density functionals. 
It consists of interacting electrons in an infinite volume 
and in the presence of a uniformly distributed background 
positive charge. Traditionally, in its paramagnetic version, 
the system is constructed by allowing the number n of 
paired electrons in a D-dimensional box of volume V to 
approach infinity with p — n/V held constant. ^''^ 

Using atomic units, the high-density (r^ — >■ 0, where 
Tg is the Seitz radius) expansion of the reduced energy 
(i.e. energy per electron) of ZJ-jellium is 

e{D) = eT{D)+e^{D) + ec{D), (1) 

where ex and ex are kinetic^ and exchang^^ energies 

£t[D) = ^ — , ex[D) = , (2) 

and £c is the correlation energy. After many important 
contributions,'21t26| jg known that, for 2- and 3-jellium, 
the correlation energy takes the form 



£c(-D) =^K(i?) + &j(i5)lnr,]^ 



(3) 



The constant term in Q is usually decomposed as 

ao{D)=ao,i{D)+ao,K{D), (4) 

where Oq^j is the direct ("ring-diagram") contribution, 
and ao,K is the second-order exchange part. The first few 
aj and hj and are known analytically or numerically for 
the important D = 2 and D — 3 cases (see Table |l| . 

In this Article, we introduce an alternative paradigm, 
in which the electrons are confined to a Z?-sphere, that is, 
the surface of a (-D -I- 1) -dimensional ball. These systems 



possess uniform densities, even for finite n, and because all 
points on a D-sphere are equivalent, their mathematical 
analysis is relatively straightforward.^^ 

Electronic properties of the UEG on a 2-sphere have 
been previously studied in modeling multielectron bubbles 
in liquid helium (see Ref. ,32J , and similarities between 
this system and 2-jellium have been noticed by Longe 
and Bose.l^ However, the UEG on a 3-sphere has not 
been considered before, and this Article presents the first 
study of correlation effects in a spherically-confined three- 
dimensional UEG. 



II. HARTREE-FOCK ENERGIES 
A. Exact results for finite n 

The orbitals for an electron on a D-sphcre of radius R 
are the normalized hyperspherical harmonics Yg^, where 
£ is the principal quantum number and /i is a composite 
index of the remaining quantum numbers.'^^^Il confine 
our attention to systems in which every orbital with 
£ = 0,1, . . . , L is occupied by two electrons, thus yielding 
an electron density that is uniform on the sphere (see 
Eq. ^ below) . The resulting model is defined completely 
by the three parameters D, L and R. 

The volume of a I?-sphere is 

where F is the gamma function,'^ the number of orbitals 
with quantum number £ is 

{2e + D-l)T{£ + D-l) 

r(D)r(^ + 1) ' 

and each of these has energy 

£{£ + D- 1) 



rii 



(6) 



2i?2 



(7) 



Because the total number of electrons is 
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it follows that the uniform electron density is 
n _ r(D/2 + l) 

{2L + D)r{L + D) 1 
r{D/2 + l)r(L + 1) (47ri?2)i)/2 ■' 



P 



(9) 



and the Seitz radius is 



T{L + l) 



-,1/D 



with 



_iL + D/2)T{L + D) 
=2i-i/^r2/^(i^/2 + l 



(10) 



(11) 



Using the hyperspherical harmonic addition theorem, '^^ 
one finds that the one-particle density matrix is 

L ni 

_ T {D/2 + 1)T{L + l)piD/2,D/2-i) 



T{L + D/2 + l) 



(cos 6*), 
(12) 



where P^'^'' is a Lth degree Jacobi polynomialP^ The 
angle 9 is that subtended by the electrons at the origin and 
is related to the interelectronic distance by the relatiorP^ 



ri2 



In -ral - 2i?sin(6l/2). 



(13) 



The density matrix decays rapidly with interelectronic 
separation when L is la rge (F ig, [l]), illustrating the "short- 
sightedness" of matter!3SESl 

Many properties of the UEG on a D-sphere can be found 
from Eqs. (|6| - (12). Its kinetic energy, for example, is 



£T(Ai) = -E 



D L{L + D) 
2{D + 2) ^2 ■■ 



(14) 



and it can be shown that its exchange energy is 



ex{D,L) = - 



pi(ni,f22)' 



"in J J ri2 

n r(^) DT{D-i) r(£ + f; 

2R^T{^) 2L + D r(L + Z3-l) 

-L, L + D, 
1 

2 ' 



D-1 _ 1 

-L-h, L + D-i 4^ 



(15) 



where 4-^3 is a generalized hypergeometric function.l^Sl 

B. The thermodynamic limit 

The above expressions are exact for all L but, in the 
thermodynamic limit (n, L — > cx)), each simplifies signifi- 



TABLE I. High-density coefficients for D-jellium and the UEG 
on a D-sphere. j3 and ( are the Dirichlet beta and Riemann 
zeta functions .ISSI 



Coefficient Term 



D = 2 



a-2 
a-i 
bo 

ao 

bi 



' s 

^0 



, In r. 



ao,j 



1/2 

_ iV2 
3n 



ln2 - 1 

m - ^m) 

-^/2(f -1) 



^(9-/4)^/^ 
-£(9vr/4)^/^ 
(1 -In2)/7r2 

-0.071099 

+0.009229 



TABLE II. Numerical values of ao,j(2, L), ao,K(2, L), co,j(3, L) 
and ao,K(3,L) for various L. 



L 




UEG on a 


2-sphere 




UEG on a 


3-sphere 


n 


ao,j(2,L) 


ao,K(2, L) 


n 


co,j(3, L) 


ao,K(3, L) 





2 


"0.2274 


+0.1137 


2 


-0.0476 


+0.0238 


1 


8 


-0.2534 


+0.1111 


10 


-0.0717 


+0.0231 


2 


18 


-0.2677 


+0.1118 


28 


-0.0897 


+0.0231 


3 


32 


-0.2762 


+0.1124 


60 


-0.1038 


+0.0233 


4 


50 


-0.2816 


+0.1128 


110 


-0.1154 


+0.0234 


00 


00 


-0.3069 


+0.1144 


00 


—00 


+0.0242 



cantly. For example, 



r(i:i + i) 

2 



L 



D 



T{D/2 + l) (47ri?2)D/2^ 
R 



Pi ^ P 



r{D/2 + l)JD/2{L0) 

{Le/2)D/2 



(16) 
(17) 
(18) 
(19) 



where J„ is the nth-order Bessel function.l^ We note that 
( 19 1 reduces to the usual density matrices in 2-jelliurrP^ 
and 3-jellium.ISThe kinetic and exchange energies become 



eT(^) = 
sx{D) = 



D 

'2{D + 2)^' 
2D an 



(20) 
(21) 



Equations (20) and ( |2T| ) yield the two terms in ([2]), and 
are identical to the Z?-jellium expressions. Particular 
cases are given m Table |I] These results were originally 
discovered by Glasser and Boersma,'^ and IwamotcP^ for 
-D-jellium, but our derivation for the UEG on a Z3-sphere 
is more compact than theirs. 
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FIG. 1. The one-particle density matrix for the UEG on a 
3-sphere. Plots for L = 5 (blue), L = 10 (red), L = 20 (yellow) 
and L = 40 (green). 



ao.K (D, L) 




(L+1)- 



FIG. 2. ""■'^7',' as a function of (L + 1)^^ for D = 2 (blue 
dots) and D = 3 (red squares) . 



III. CORRELATION ENERGIES 

A. Exact results for finite n 

We now turn our attention to the study of the cor- 
relation energy of the spherically-confined UEGs. By 
applying perturbation theory to UEG on a 2-sphere, we 
find that the reduced energy coefficient corresponding to 
the lowest-order ring-diagram contribution is 



occ virt 



ao 



71 ^ Ki -\- t^a 

I, J a,b 



coefficient in S0(4) is much simpler than in S0(3^ and 
the energy coefficient from the lowest-order ring-diagram 
is 



^ L oo 

co,j(3,L) = -^ E 



(z + l)(j + l)(a + I)(6+l) 



i,j=Oa,b=L+l I J a 



E 



(2/7r)2 



(^-t-i)2(£+|)2^ 



(23) 



where the sum over i respects the same restrictions as in 
the 3j symbols in ( 22 ) . 



The second-order exchange part for the UEG on a 
2-sphere is 

{ij\ab){ba\ij) 

ao,K 2, L)^- y^y^ — , 



E E 



{2t + l){2j + l){2a + l){2b + I) 



» j=0 a,b=L+l Ka + Kh-Ki- 

" f-; looo/looo/ 



,0 / lo oi Ij £2 M ' 



(24) 



where the curly brackets denote 6j symbols, '^^ and for the 
UEG on a 3-sphere, we found 



ao 



k(3,L) = -E E (^ + ^)(-^' + ^)(" + ^)(^+^) 

'^*j=Oa,f,=L+l Ha+ Hb- 

(2/^)(4 + l) (2/^)(^2 + l) \i h a\ 



E 



(^1 + + i) (^2 + i)(^2 + i) \j i2 bj' 

(25) 

where we have used the S0(4) version of the 6j symbols.!^ 
Numerical values for finite L are given in Table [llj 

B. The thermodynamic limit 



In the thermodynamic limit {n — >■ 00), Eq. (22) be- 
comes 



1 ^ j2 (2i + l)(2j + l)(2a + l)(2b+l) 



^ 2e + i 



1 ( i e a\ j e b 



,0 00/ \o 0, 



(22) 



ao,j(2) = - 



where 



00 



d£ , , /■^+' 
— / idt jd] 

-Iri -In Jmax(r,|£-i|) 



a da 

Jo Jo 

X 1'''' bdb „ -^'jl-^'^f'" (26) 



max(l,|^-j|) 



fj2 _|_ ^2 _ j2 _ j2 



where {ij\ab) are two-electron integrals and the brackets 
are 3j symbols!^ For the UEG on a 3-sphere, the coupling 



TT ^{i + £ + a){i + 1- a){£ + a - i){a + i - £) 

(27) 
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comes from the asymptotic expansion of the 3j symbolP 



Defining a 



2iicos9i and 6^ 



2ji cos 62 to transform the a and b integrals into 6*1 and 
02 integrals, and then renaming i, j and £ as pi, p2 and 
q, we obtain 



1 

47r3 



7^ 



|pil<i 
|pi+q|>i 



dpidp2 



IP2K1 q2 

|P2+q|>l 



(28) 



which is precisely the lowest-order ring-diagram contri- 
bution of 2-jelliumfll' This integral can be solved^'' to 



yield 



ao,j(2) =ln2-l. 



(29) 



One also finds that the higher-order ring-diagram contri- 
butions are identical in 2-jellium and in the UEG on a 
2-sphere and the resummation technique^^ yields 



&o(2) = 0, 



6i(2) 



10 

3^ 



1 



(30) 



For large L, the sums in Eq. (23) can be replaced by 



integrals and one finds and the leading order of ( 23 1 is 



co,j(3,i) 



d£ 

l/L Jo 
1 ri+e 



dt 



In 2 



ai 



(■''"-'"^'dadi 



(31) 



In- 



L 



1 -ln2 



Inr, + 0(r" 



ft follows that 



60 (3) = 



l-ln2 

9 ' 



(32) 



and thus the logarithmic divergence of the correlation 
energy in the UEG on a 3-sphere is exactly the same as in 
3-jellium. One notes that the result (31 1 can be derived 
for any value of the radius R. The latter divergence, 
contrary to some claims in the literature, does not result 
from the long-range part of the Coulomb operator but 
from its short-range part.*'"' The observation of the 
same divergence in the UEG on a 3-sphere — where 
the interelectronic distance can never exceed 2R — also 
demonstrates this. 

Proceeding similarly to the D = 2 case, it can be shown 



that, in the thermodynamic limit, (23) becomes identical 
to the expression of the second-order ring-diagram in 
3-jellium: 



co,j(3) = 



167r5 



IpiKi 
|pi-i-q|>i 



dq 



dpidp2 



IP2KI q2 
|P2+q|>l 



- q ■ (pi +P2) 



(33) 



where the excitation vector q has the domain < q < 
00!^ Moreover, the higher-order ring diagram contribu- 
tions are also identical in 3-jellium and in the U EG on a 
3-sphcre. Using the resummation techniquep^ it follows 
that-^ 



ao,,i(3) = -0.071099. 



(34) 





-.m - - 


fO. 114357, 


In 2 






IT " 




+0.024179, 



For D — 2 and 3, we have not been able to prove the 
equivalence of the second-order exchange contributions 
in D-jellium and in the UEG on a Z?-sphere, but the 
numerical results in Table [ll] and Fig. [2] suggest that, in 
the thermodynamic limit, ao,K(2) ~ +0.11 and ao.K(3) ~ 
+0.024, which may be compared with the known 2-jellium 
and 3-jellium values;^ 



ao,K(3) 



where /? and ^ are the Dirichlet beta and Riemann zeta 
functions .I^Sl 



IV. DISCUSSION 

Uniform electron gases on a D-sphere are an attractive 
generalization of D-jellium and, as we have shown, one can 
derive compact expressions for the first few terms of the 
high-density energy expansions for both finite and infinite 
systems. Although UEGs on a _D-sphere are physically 
different from D-jellium, we have shown that, in the 
thermodynamic limit, the first few terms of their high- 
density energy expansions are identical and we conjecture 
that the high-density expansions are identical to all orders. 

Recent calculations on the Thomson problem suggest 
that the leading term of the low-density (large-r^) energy 
expansions in 2-jellium and in the UEG on a 2-sphere are 
also identical.'^ Moreover, because the Thomson problem 
is trivial for Z? = 1, it is actually possible to show the 
strict equality of the leading term of the low-density 
energy expansions in 1-jellium and in the UEG on a ring 
(l-sphere).E2l 

Although it is pleasing to know that the spherical and 
conventional gases become equivalent in the thermody- 
namic limit, we believe that it is even more important 
to recognize that they are not equivalent for finite n. 
Equations (14) and (15) predict significantly different ki- 



netic and exchange energies from (20) and (21 1 when n 
is small. Moreover, combining the information from the 
high- and low-density regimes, one can easily construct 
local-density approximation-type correlation function als 
for finite systems using interpolation functions.'^^'i^^ We 
believe that the UEG on a D-sphere will be useful in 
the future development of correlation functionals within 
density-functional theory.*^ 
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